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ABSTRACT 


To use H gas as a common fuel it needs to be in high pressures or cryogenic temperatures to 
have reasonable density. But, if we have adsorbent materials with high volumetric capacities to 
store hydrogen at ambient temperature and low pressure without any compressing it is 
worthwhile to use hydrogen as clean and reversible fuel. Here, we want to report the 
adsorption and decomposition properties of aluminum and gallium-doped graphene at 
ambient temperature. We studied the adsorption of H20 molecule on pure and doped 
graphene via density functional theory. So, possible interactions between the HO molecule 
from three sides and pure and aluminum and gallium-doped were examined. After adsorption, 
decomposition of the HzO molecule has been studied and so on, for receive a reaction 
pathway, possible intermediates and transition states has been studied. To continue the 
density of states, interaction energies and thermodynamic parameters have been calculated. 
The results showed that the adsorbed water on aluminum and/or gallium-doped graphene 
decompose to OH and H and then adsorb on the surface again at ambient temperature and 
this process was thermodynamically favorable. 
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INTRODUCTION 


One of the useful ways for energy storing is the chemical bonds. The advantage of this way was known for decades 
as the combustion of fossil fuels and breaking C-H bonds. However CO,, NOx, and SO, emissions from fossil fuels 
cause to necessitate a transition to non-carbon energy fuels [1]. Hydrogen with the energy density about of three 
times more than gasoline and water as it's the only waste product would help us in the transition to a non- 
emission energy system [2]. Indeed, it could encourage application of renewable energy and help transform to use 
non-carbon fuels. Today using hydrogen for transportation as hydrogen fuel-cell are operating on the road and it 
is already realized with over 120 hydrogen stations in California and Germany [3]. So, this is encouraging to use 


+Corresponding author: s.fshari@du.ac.ir 
1 of 13 
This is an open access article published under the CC BY 4 DEED license 


Materials Chemistry Horizons 2024, 3, 1056 2of 13 


hydrogen as a transportation fuel, but this replacement is limited by the problems of hydrogen storage, Finite 
space of the vehicle and especially the requirement to systems that enhance driving ranges need high density and 
volume of energy, respectively. On the other hand, storing the hydrogen as a liquid needs a very low temperature 
with densities of 70 g/L and as a compressed gas needs a very high pressures (more than 700 bar) with densities of 
39 g/L [4]. Refueling by both ways is expensive and also developing of storage systems for operating under these 
excessive conditions is difficult. Furthermore, the storage of hydrogen at that pressure needs significant energy for 
compression. Also, storing at very low temperature increases storage capacities but seriously reduces the system 
volumetric density and causes costs for operating the thermal management. So, we need high-density hydrogen 
storage with moderate conditions and no cryogenic or highly compressed hydrogen storage [5, 6]. Adsorbent- 
based systems could be a suitable alternative for this aim. Hence, we can use materials capability for storing at 
ambient temperature. In the employing tanks for the hydrogen gasses storing the operating at low pressures, 
being light, less expensive, and more conformable parameters are serious. So, by these or similar capabilities, the 
fuel cell using in transportation will be provided as well as other usage [7]. 

In recent years, many scientists interested into design and synthesis of new materials with mentioned capabilities 
for hydrogen storage [8-12]. Nanomaterials have been investigated as promising adsorbents for this. In particular 
the high surface-to-volume ratio and chemical versatility of these materials cause to increase hydrogen uptake. It is 
important to extend active, stable non-toxic and inexpensive surfaces, too. While, using of the transition metal- 
based adsorbent are limited due to cost and toxicity, there are many studies that report, graphene displays an 
acceptable surface activity [13-15] and displays a good performance as a catalyst for dehydrogenation reaction 
[16-18]. Water adsorption on nanostrucures like graphene is interested in for many researchers. Amanda J. Carr et 
al., studied on the adsorption of the water on graphene and graphene oxide experimentally [19]. Ikutaro Hamada 
studied about the adsorption of water on graphene by the van der Waals density functional study and compared 
the results with the quantum Monte Carlo approaches [20]. In addition, many studies have approved the 
modifying of the physical and chemical properties of graphene by substituting the carbon by the heteroatom. It is 
because of the changes in electronic structure [21-23]. For example, the surface reactivity of graphene and other 
nanostructures could be improved by the aluminum and gallium doping [24-26]. Also, Yong-an Ly, et al. showed 
that the aluminum or gallium doping enhanced the graphene activities for adsorption and dissociation of N20 
[271. Qilin Yi, et al. showed the aluminum or gallium doping improved the graphene surface adsorption properties 
for NF3 and NO molecules [28]. In addition to the storing problem, we need to have an iminium source of 
hydrogen to produce and then storage. Hydrogen molecule can be produced in different ways. Decomposition of 
hydrogen-rich material is one of the most attractive ways for hydrogen-production. One of the important 
hydrogen sources for this way is hydrocarbons. Water is the other hydrogen-rich, Water is a safe material respect 
the hydrocarbon and is the most material on the earth. For this, water decomposition was interested for many 
scientists [29, 30]. So, by this way, hydrogen atoms could be produced from the water adsorbed and stored on the 
surface and then desorb to produce the hydrogen gasses for use. 

In the preset work, the adsorption and decomposition of the H20 molecules on the Al and/or Ga-doped graphene 
sheet have been studied in room temperature and using density functional theory (DFT) calculation. The main aim 
of this study was the investigation of the HO molecule's decomposition mechanism on the graphene surface. 
Also, the possible reaction pathways and intermediates have been studied for this process. The performance of 
heteroatom doping on the potentials of graphene for this process and use as hydrogen storage has been studied, 
too. 


COMPUTATIONAL METHOD 


In this work, the investigation on the HO molecules to adsorption and subsequently dehydrogenation on the 
graphene sheet were performed by using density functional theory (DFT) and time depending density functional 
theory (TD-DFT) and M062X/6-31g(d) method/basis sets. For this study, we have considered the armchair graphene 
sheet consisting of 25 carbon honeycomb that the end carbon atoms were bonded and saturated by hydrogen 
atoms. We have studied non-pure graphene, too. So, we have considered the Al-doped graphene as configuration 
"#,, Ga-doped graphene as configuration "8%, 2Al-doped graphene as configuration "C,, 2Ga-doped graphene as 
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configuration "D” and, Ga, Al doped graphene as configuration "E”. We have used the QST3 method to calculate 
the transition state (TS). To confirm the stationary structure with one imaginary frequency, we have done the 
frequency calculations. All the calculations have been done using the GAMESS program suite [31]. To illustrate the 
configurations, we used GaussView 5.0, 

We have used the following equation to calculate adsorption energy (Eas) of adsorbed HO molecule(s) on Ga 
and/or Al-doped graphene sheet, 


Bay 


E(surface+water) — (Esurface + Ewater) a) 


where Eur fax iS the calculated energy of the Ga and/or Al-doped graphene sheet, Eater is the calculated energy 
of one HO molecule and Esur face+water iS the calculated energy of one HO molecule adsorbed on Ga and/or Al- 
doped graphene sheet. 

Also, we have studied the density of states, HOMO, and LUMO energy levels for the configurations before and after 
possible interactions, For this, we used the GaussSum program [32]. 


RESULTS AND DISCUSSION 


In this work, we have studied the adsorption and decomposition of the HO molecules on Ga and/or Al-doped 
graphene sheet. All structures of the configurations were optimized, energetically at ambient temperature. Also, 
the surface activity was studied for the doped graphene for Hj0 molecules adsorption and then the 
decomposition of two HO molecules. The size of the considered configurations [33] and the level of theory [34] 
employed for this study were validated by previous works. To continue we will discuss the process of adsorption 
and dehydrogenation of 0 molecules on the mentioned surfaces, 


Adsorption 


Figure 1 illustrated the structure from top and side view for Al-doped graphene (configuration ‘A")) when M = Al or 
Ga-doped graphene (configuration "B")) when M = Ga before and after adsorption of one H,O molecule. Also, the 
optimized structure of 2Aldoped graphene (configuration °C") when M1, M2 = Al, 2Gadoped graphene 
{configuration "D")) when M1, M2 = Ga or Ga, Al-doped (configuration "E")) when M1 = Ga and M2 = Al, before and 
after adsorption of two water molecules from top and side view have been shown in this figure, 

As is clear in Figure 1, the optimized structures of mentioned configurations after adsorption of the water 


molecules changed and be out of the plane, while, they were two dimensions sheetlike graphene before 
interaction, For study about the adsorption of HO molecules on graphene, the HO molecule is set upright on 
graphene from H atom, from O atom, and one side of the HO molecule. Adsorption has been studied on the 
carbon atom and one of the honeycombs, For prevent of the edge effects, the center of the studied graphene was 
selected. There was no any interaction between the HO molecule and the mentioned studied sites of graphene. 
On the other word, in this work, pure graphene has not adsorbed any water molecule. Then, one of the carbon 
atoms of the graphene has been substituted by one Al atom (configuration A).. For this configuration, adsorption 
process has been studied on the doped atom site and the HO molecule from three mentioned sides. It was seen 
the only adsorption from the O atom of the HO molecule, and there was no adsorption from two other sides on 
the Al atom for configuration A. In continue configuration 8, that one carbon atom of graphene has been 
substituted by Ga, was designed. For this configuration, the adsorption process has been done as the same as 
configuration A, and the HO molecule only was adsorbed from O atom on Ga atom, too. As can be seen from 
Table 1, the scale of adsorption energy for the studied molecules is in the chemical adsorption range. 

The density of states (D0S) was calculated for graphene and configurations A and B before and after one HzO 
molecule adsorption. The DOS diagrams for the graphene and mentioned configurations before and after 
adsorption have been shown in Figure 2. The energies of the HOMO, LUMO levels, and the difference between 
them (Ey) for these configurations have been shown in Table 1. As is clear from Table 1 and Figure 2, by the 
adsorption, the HOMO and LUMO energy levels were shifted and also the Ey increased. The E, for our studied 
graphene was 2.00 eV, which is in good agreement with other previous studies [35-38]. By the substitution of one C 
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atom by Al or Ga atom the Ey decreases to 1.77 and 1.78 eV, respectively. It means that by doping the conductivity 
of graphene increases. After adsorption of one H,0 molecule on both configurations A and B, the Ey were 
changed to 2.21 eV. That change was about 0.43 and 0.42 eV for configuration A and B, respectively. It shows that 
both the Al-doped and Ga-doped graphene are sensitive to water adsorption, 
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Figure 1. The optimized structure of a) Aldoped graphene (M=A) or Ga-doped graphene (M=Ga) on top (left) and side (right) view, 
by Al or Ga-doped with an adsorbed water molecule on top (left) and side (right) view, c} 2Aldoped graphene (M1, M2=A\), 2Ga- 
doped graphene (M1, M2=Ga) or Ga,Abdoped (M1=Ga and M2=A\) on top (left) and side (right) view, d) 2Al, or 2Ga, or Ga,Ardoped 
graphene with two adsorbed water molecules on top (left) and side (right) view. 

To continue, in the graphene sheet, two carbon atoms were substituted by two Al atoms (configuration C).. For this 
configuration, adsorption process was studied for the mentioned doped sites and it was seen as the only 
adsorption from © atom of the water molecule. Then configuration D is designed which two previous carbon 
atoms were substituted by two Ga atoms. For this configuration, the adsorption process has been done as the 
same as configuration C. The last configuration was € that two carbon atom were substituted by one Ga atom and 
one Al atom. Also, for this configuration, the HO molecules only were adsorbed from the O atom on Ga and Al 
atoms. 

The DOS was calculated for configurations C, D, and € before and after adsorption of two HO molecules, too. The 
DOS diagram for the graphene and three mentioned configurations before and after the adsorption of two HO 
molecules have been shown in Figure 3. The energies of the HOMO, LUMO levels, and E, have been shown in Table 
1, too. As the Table 1 and Figure 3, the HOMO and LUMO energy levels were shifted and also the E, increased for 
configurations C, D, and E after the adsorption. In Figure 3 (a),, the DOS of graphene have plotted beside the DOS 
of configurations C, D, and E for comparison, By the substitution of two C atoms by Ga and/or Al atom to design 
the configurations C, D, and E the Eg decreases to 1.89, 1.95, and 1.70 ev, respectively. So, for these configurations, 
the conductivity is more than graphene, too. The E, for the configurations C, D, and E was changed to 1.82, 2.02, 
and 2.01 e¥, respectively by the adsorption of two water molecules. These changes were decreasing for 
configuration C by about 0.07 eV and increased by about 0.07 and 0.31 eV for configurations D and E, respectively. 
It shows that configuration E is more sensitive than configurations C and D to water adsorption. Also, the amount 
of charge transferring (qcr) from water molecule to the Al and/or Ga-doped graphene after adsorption process was 
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listed in Table 1, As can be seen from Table 1, the amounts of qCT for all interactions between water and Al and/or 
Ga-doped graphene are positive, 


Table 1. The adsorption energy, energies of the HOMO, LUMO levels, &, and qcy for the graphene, configurations A, B, C, D, and 
Ebefore and after adsorption of ane or two water molecule (s). 


Configuration Egg (kcal/mol) HOMO (ev) LuMo (ev) Fev) acre) 

Graphene 3 479 2.79 2.00 2 

A : “4.93 a6 wz : 
A+H,O na 4.80 2.59 221 0.16 

B é 4.95 317 1.78 : 
B+H,0 77.35 4.80 2.59 221 0.15 

c : 4.86 297 1.89 : 
c+2H,0 129.01 4.23 241 1.82 0.16 

> : 4.89 2.94 1.95 : 
D+2H,0 4141.52 4.38 2.36 2.02 031 

E 4 ABA 314 1.70 : 
E+2H,0 4137.53 4.36 2.35 201 0.32 

Decomposition 


Decomposition of One H20 Molecule 


After adsorption of the HO molecule on the configuration A (or B), the decomposition process starts, For 
decomposition of H20 molecule it planned to be decomposed two O-H bonds. First step is involving the stretching 
and then breaking of the two O-H bonds. As configurations A and B are the same for the decomposition steps, we 
only showed the decomposition process of one HO molecule on configuration Ain Figure 4, for briefness. But, for 
both configurations the thermodynamic parameters (reaction energy (AE), change of Gibbs free energy (AGzog) and 
change of enthalpy (AHy93)), imaginary frequencies (v), and corresponding activation energies (E,q) have been 
summarized in Table 2. 


Table 2. Calculated parameters for every possible steps of decomposition of one water molecule an configurations A and 6. 
(activation energy (Eq), imaginary frequency (v), reaction energy AE, change of Gibbs free energy (AGy9g), and change of enthalpy 


(Hoa). 
Tee Fact v aE G58 Hao 
(kcal/mol) (kcal/mol) (kcal/mol) {kealmol) 
Configuration A 
P12 145 1350) 128 129 134 
P24P4 22.4 2951 20.9 a7 818 


Configuration B 


PIoP2 19.7 13761 Ad 41 45 
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Figure 2. 005 diagrams for a) graphene, configuration A and configuration B, b) configuration A before and after adsorption of 
‘one water molecule, and ¢) configuration B before and after adsorption of one water molecule. 


Materials Chemistry Horizons 2024, 3, 1056 7 of 13 


Figure 3. DOS diagrams for a) graphene, configuration C, configuration O and configuration €, b) configuration € before and after 
‘adsorption of two water molecule, ¢) configuration D before and after adsorption of two water molecule, and d) configuration € 
before and after adsorption of two water molecule. 

Also, Figure 5 have shown the decomposition reaction pathways for a HO molecule on both configurations A and 
B. In these ways, the first intermediate (P1) was taken as one adsorbed H20 molecule on the configuration A and 
configuration B (see Figure 5). Breaking of the O-H bond in first step of the reaction pathway starts with the 
elongation of it in the adsorbed HO molecule. In the first step, for the mentioned configurations, the energy 
barrier was 14.5 and 19.7 kcal/mol, respectively. For the configuration A this is smaller than the corresponding 
activation energy over Co (110) surface (ca. 19.14 kcal/mol) [39]. It shows that these reactions is thermodynamically 
favored and so takes place at ambient temperature. As can be seen from Figure 4, at the TS1 (first transition state) 
the 0-H2 bond elongated from 0.97 to 1.21 A (configuration A),, and 1.25 A (configuration B).. Also, the O—Al and 
(O—Ga distance’s was shortened from 1.94 to 1.85 and from 2.00 to 1,92 A, respectively. By elongation the O-H2 
bond the second intermediate (P2) has been created. The lengths of newly formed C1-H2 bond in the second 
intermediate are 1.11 A, for both of the mentioned configurations, The energy of the P2 for the configurations A 
(-12.8 kcal/mol) and B (-4.1 kcal/mol) is lower than the P1 
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Figure 4. The stationary points structures for dehydrogenation af H20 molecule on Al-doped graphene for configuration A after 
‘energeticaly optimization, 


Initial adsorption energy -71.41 (configuration A) and -77.35 kcal/mol (configuration B) and dehydrogenation 
barrier for configurations A and B are 14.5 and 19.7 kcal/mol, respectively. As can be seen, in this way the initial 
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adsorption energy is large enough to the subsequent dehydrogenations. 
In the second step, by breaking another hydrogen atom the dehydrogenation has been studied. For the P2, O-H1 
bond could break while a large activation barrier of about 82.4 kcal/mol for configuration A (Figure 5), and for 
configuration 8 this step did not occur: On the configuration A, the TS1 structure is verified by an imaginary 
frequency of the 0-H2 stretching vibration mode (1350i cm”), On the configuration B, the imaginary frequency for 
the same stretching vibration mode is 12761 cm”, Also, for the configuration A, the TS3 structure that is 
correspond to the second H atom on the HO molecule for dehydrogenation (O-H1) is verified by an imaginary 
frequency (295i em~), 

‘As can be seen from Tables 1 and 2, the first step with AHyog = -13.4, AGyog = -12.9 kcal/mol, and reaction energy of 
-71.41 kcal/mol respect to the initial state for configuration A, and with AHy9g = ~4.5, AGagg = ~4.1 kcal/mol, and 
reaction energy of -77.35 kcal/mol respect to the initial state for configuration 8 is exothermic and 
thermodynamically favored. While, the second step with AHog = +81.8 and AGz9q = +81.7 kcal/mol is endothermic 
and is not thermodynamically favored. It shows that for the adsorbed HO molecule on both configurations A and 
B, the second step of dehydrogenation is not occurred easily at ambient temperature. 


Decomposition of Two H20 Molecules 


To continue, the decomposition of two adsorbed HO molecules on the configurations C, D, and E have been 
studied. The same dehydrogenation steps have been studied for these configurations and so, only the reaction 
pathway for decomposition of the two adsorbed HO molecules on configuration D and the optimized structure of 
stationary points have been shown in Figure 5, for briefness. So, after adsorption of two H20 molecules, as 
mentioned, onto the Al atoms (configuration C), or Ga atoms (configuration D),, or Ga & Al atoms (configuration 
), the decomposition would be start. The thermodynamic parameters (change of Gibbs free energy (AG), reaction 
energy (AE), and change of enthalpy (AH), imaginary frequencies (v), and corresponding activation energies (Esq) 
have been summarized in Table 3. 

In Figure 5, the first intermediate (P1) has been taken as configuration D with two adsorbed HO molecules. By 
studding on the first HO molecule with elongation of the O1-H1 bond, the reaction pathway of dehydrogenation 
begins. As can be seen in Table 3, the activation energy barrier for this step is 21.8 kcal/mol. This energy barrier is 
13.5 kcal/mol for configure C and 16.3 kcal/mol for configurations E. This shows that these reactions would occur 
rapidly at ambient temperature. The O1-H1 bond elongated from 0.97 to 1.23 A, and the O—Ga distance was 
shortened from 1.99 to 1.94 A at the first transition state (TS1) on the configuration D (see Figure 6). By elongation 
the O1-H1 bond the second intermediate (P2) has been created. The length of newly formed C1-H1 bond in the P2, 
is 1.11 A, and the Gai—O1 distance's stays at 1.84 A, The energy of the P2 for the configuration D is -7.1 kcal/mol 
and is lower than the P1 and also, the initial adsorption energy for this configuration is about -141.52 kcal/mol 
and dehydrogenation barrier energy is 21.8 kcal/mol, It is clear that, in this way the initial adsorption energy is 
large enough to the subsequent dehydrogenations. 

In the next step, breaking another hydrogen atom for the dehydrogenation process has been studied. On 
intermediate P2 configuration the activation barrier energy for 02-H2 bond splitting was about 27.4 kcal/mol (see 
Figs. and 6). As the previous step and the first adsorption process are exothermic, this energy barrier could partly 
offset. On the surface of configuration D the TS1 structure is verified by the imaginary frequency that related to the 
stretching vibration mode of O1-H1 bond and is 14981 cm”, For the TS2 configuration the bond elongation 
change for 02-H2 is from 0.97 to 1.32 A, and the O—Ga2 was shortened from 2.0 to 1.90 A. Also, the TS2 structure 
is verified by the imaginary frequency that related to the stretching vibration mode of 02-H2 bond and is 11541 
cm”, In the next step, it needs to be study about intermediate 3 that is obtained by elongating and breaking the 
02-H2 bond. The length of newly formed C3-H2 bond in the P3 configuration is 1.11 A, and the distance of Ga2—- 
02s shortened from 2.0 to 1.81 A. 
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Table 3. Calculated parameters for every possible steps of decomposition of two water molecule on configurations C, D and E 
(activation energy (Eac), imaginary frequency (v), reaction energy AE, change of Gibbs free energy (AGz9q), and change of enthalpy 


(ahiz09) 
Reaction oa S oe ed mae 
(eatvmo) em") (kealmol) (kcal/mol) (kcal/mol) 
Configuration ¢ 
P12 135 13581 312 288 327 
P2473 35.1 328 40 28 53 
p34 174 5961 oa.e 70.0 683 
Paps 655 3761 607 60.2 605 
configuration D 
P12 218 14981 a ana 22 
P24P3 274 114i 37 30 40 
P34 656 16991 530 589 S73 
Paps a8 nei 620 378 576 
Configuration E 
P12 163 13131 55 446 487 
P24P3 268 11631 49 28 34 
P3.P4 n3 1677 615 629 619 
Paps 234 10591 556 38.7 378 


Figure 5. One H20 molecule decomposition reaction pathways on configuration A and 8, and two H20 molecules decomposition 
reaction pathways on configuration C, D, and E. 

By the elongation of the O1-H3 bond, study about the decomposition reaction pathway of adsorbed water 
molecule continued. In the TS3 configuration the O1-H3 band was elongated to 1.37 A, and the 01—Ga1 distance 
was shortened to 1.87 A. An imaginary frequency of 16991 cm~" in the TS3 is related to the mode of O1-H3 bond. 
Then, for obtained the intermediate P4 the O1-H3 bond breaks. The length of C2-H3 bond in the P4 is about 1.10 
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A, and the Ga1—01 distance stays at 1.87 A. The reaction energy of this step is about +53.0 kcal/mol and it needs 
457.1 and AGaog = 458.9 
kcal/mol. So, this reaction would not arise at ambient temperature because this is not thermodynamically favored. 

In continue of decomposition process, the 02-H4 bond at the TS4 configuration was elongated to 1.37 A, and the 
(02—Ga2 distance was shortened to 1.85 A and an imaginary frequency of 718 cm’ in the TS4 is related to 02-H4 
bond elongation. By breaking, the 02-H4 bond, intermediate PS is obtained. The length of C4-H4 bond in the PS 


to pass the activation energy barrier about 65.6 kcal/mol, relative to P3, with AHyog 


configuration is 1.11 A, and the Ga2—02 distance stayed at 1.87 A. The calculated activation energy and reaction 
energy for this step is 91.8 and +62.0 kcal/mol, respectively. Also, as can be seen from Table 3 AH and AG for this 
step is +57.6 and +57.8 kcal/mol, respectively. These results mean that this reaction would not arise at ambient 
temperature because of the thermodynamic limits. So, for this configuration, the second H atoms of two water 
molecules (H3 and H4) bonded to 01 and 02 atoms, respectively, instead to bond to the carbon atoms of studied 
surface. This is because the aromaticity of the sheet. The O atom creates double bond with Ga atom of 
configuration D and it cause to disturb the aromaticity of the sheet. So, it means that the reaction is not kinetically 
preferable 


oumisa 
CEH 1484 
O18. 097A 
GaLOir194A, 
Ga202: 2004 
02-2, 097A 


o1an:266 
Cinta 
O1BB:097A, 
Gab Or isa, 
a2-02 2004 
oan: 097 A 


OLEL26 A 
Ceara 
Lan. 0s7A 
Gal-O1- 1838, 
Ga202: 1.908 
OL OSTA 
can: 137A 


onm:261 
CLHILIA 
o1n-097A 
Gui-01 1834, 
GeO 181A, 
OnHOgTA 


OL 3408, 
Cra1ia 
Onis 1374 


CHB. 
GaLO1 1674 
GaLOr 183A 


OH 097A 
Ova 281A, 
cA 


OLHL:329A, 
ceria 
Oui 2174 
@w1ud 
GLO 188, 
CLONES 
Onaeis7A 
ona 3224 
CaEDTA 
CoH 126A 


Figure 6. The stationary points structures for dehydrogenation of H20 molecule on 2Gadoped graphene (configuration 0). All 
distances are in A. 

Reaction pathways for mentioned studied configurations show that the one-step dehydrogenation of the 
adsorbed H,0 molecules could be suitable. The HO molecules are decomposed to H and OH on the Al-doped, Ga- 
doped, and Al, Ga -doped graphene. In other word, after decomposition the water molecule on the mentioned 
surfaces, the hydrogen atoms is adsorbed on carbon sites on graphene sheet and OH is remained on the doped 
atoms in the studied surfaces. The H atoms could decompose and produce the Hy molecules at ambient 
temperature, while, it not possible for OH to decompose from the doped sites in this situation. It is predictable, 
because the bond energy of AI-O (119.74 kcal/mol) and Ga-O (89.38 kcal/mol) are stronger than C-H bond (80.78 
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kcal/mol) [40]. There are almost infinite sites on the studied surfaces to bond with the decomposed H atoms. So, 
we believe that the Ga and/or Al-doped graphene could be suitable catalyst for the dehydrogenation of H;0 
molecules and also convenient storage to store and then produce Hy molecules on a large scale at ambient 
temperature. 


CONCLUSION 


The pursuit of the structures capable to adsorb and storage hydrogen has already focused on storage by binding 
hydrogen on their infinite atoms at ambient temperature. We have targeted graphene as a nanostructure for its 
infinite carbon atoms to be adsorbent sites and water molecules as hydrogen-rich safe material. Pure graphene has 
not adsorbed any water molecule. So, graphene adsorption has been improved by substituting the percent of 
carbon atoms with gallium and/or aluminum atoms, Adsorption and then decomposition of water molecules have 
been done on doped graphene. While one of the water's hydrogen atoms dissociated and adsorbed on the 
graphene'’s carbon atom and the second hydrogen atom preferred to be with the oxygen atom of water bound on 
gallium or aluminum atom. This decomposition and hydrogen adsorption process was exothermic and 
thermodynamically favored, So we have an adsorbent with high volumetric capacities to store hydrogen at ambient 
temperature and low pressure without any compressing and use energy. 
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